Fast microwave beam splitters from superconducting resonators 



M. Haeberlein,i'2,[*]D. Zueco,^'^ p. Assum/ T. Weifil,^ E. Hoffmann/' ^ B. 

Peropadre,'' J.J. Garci'a-RipoU,'^ E. Solano,'''^ F. Deppe,^'^ A. Marx/ and R. Gross^-^ 

^ Walther-Metflner-Institut, Bayrische Akademie der Wissenschaften, D-85748 Garching, Germany 
^ Physik-Department, Technische Umversitat Munchen, D-85748 Garching, Germany 
^Instituto de Ciencia de Matenales de Aragon y Departamento de Fisica de la Materia Condensada, 
GSIG-Umversidad de Zaragoza, Pedro Cerbuna 12, 50009 Zaragoza, Spam 
^Fundacion ARAID, Paseo Maria AgusUn 36, 50004 Zaragoza, Spam 
^Institut Neel, CNRS, F-38042 Grenoble cedex 9 
^Instttuto de Fisica Fundamental, IFF-GSIC, Serrano 113-B, E-28006 Madrid, Spam 
"^Departamento de Quimica Fisica, Universidad del Pais Vasco UPV/EHU, Apartado 644> 48080 Bilbao, Spain 
^ IKERBASQUE, Basque Foundation for Science, Alameda Urquijo 36, 48011 Bilbao, Spam 

(Dated: February 5, 2013) 

Coupled superconducting transmission line resonators have applications in quantum information 
processing and fundamental quantum mechanics. A particular example is the realization of fast beam 
splitters, which however is hampered by two-mode squeezer terms. Here, we experimentally study 
superconducting microstrip resonators which are coupled over one third of their length. By varying 
the position of this coupling region we can tune the strength of the two-mode squeezer coupling from 
2.4% to 12.9% of the resonance frequency of 5.44 GHz. Nevertheless, the beam splitter coupling 
rate for maximally suppressed two-mode squeezing is 810 MHz, enabling the construction of a fast 
and pure beam splitter. 



Recent advances in quantum electrodynamics with su- 
perconducting circuits (circuit QED) allowed for the ex- 
perimental implementation of basic quantum computa- 
tion algorithms fP. Based on important results such as 
single photon generation ^ and multi-qubit gates O H] , 
advanced schemes for quantum error correction [5; and 
quantum feedback control were proposed. In such dig- 
ital approaches, superconducting quantum circuits sub- 
stitute classical bits and bus systems, allowing one to 
construct a general purpose quantum computation de- 
vice. However, digital quantum simulations typically re- 
quire a large number of qubits and sophisticated error 
correction schemes |7j , which is still a significant techno- 
logical challenge to date. Therefore, in the short term 
it is more promising to focus on what is called analog 
quantum computation or simulation. In this approach, 
a model quantum system is used to set up a quantum 
mechanical evolution similar to the physical system of 
interest. However, contrary to the physical system, the 
input and output channels of the model system are easily 
accessible. Superconducting quantum circuits interacting 
with quantum microwave fields represent a particularly 
attractive model system J8 . If the microwave fields are 
confined inside cavities, proposals and early experiments 
towards the simulation of manybody Hamiltonians ex- 
ist Olin]. Beyond that, recent work on systems involv- 
ing propagating quantum microwaves [HI [12] suggests 
that it is possible implement all-optical quantum simu- 
lation schemes fIBj in the microwave regime. This route 
seems particularly attractive, because superconducting 
circuits offer extraordinarily large nonlinearities jl4) and 
therefore promise deterministic gates. A qubit can, for 
example, be encoded in an entangled state of two spa- 



tially separated superconducting waveguides. In such a 
situation, linear microwave beam splitters play an im- 
portant role for the realization of single qubit rotations 
and two qubit Knill-Laflamme-Milburn gates [TSl [TC] . At 
this point, it is important to consider decoherence ef- 
fects. In order to minimize them, a beam splitter should 
be fast in the sense that its coupling rate is a signif- 
icant fraction of the frequency of the propagating mi- 
crowaves. In such an ultrastrong coupling scenario, it is 
well-known [T31[T7] that nonlinear effects arise for dipolar 
coupling. Hence, these nonlinear effects must be taken 
care of in order to ensure a pure beam splitter function- 
ality. In this work, we first develop a theoretical model 
for fast and pure microwave beam splitters based on two 
frequency-degenerate coupled superconducting transmis- 
sion line resonators with low external quality factors. We 
confirm this model with proof-of-principle experiments 
using microstrip resonators with a resonance frequency 
of ujo/2'iT = 5.44 GHz and medium quality factors ranging 
between 150 and 600. Notably, we reach a beam splitter 
coupling strength of above 800 MHz while suppressing 
the nonlinear coupling by a factor of six by exploiting 
the 90° phase shift between the inductive and the capac- 
itive coupling channel. This allows for many operations 
within decoherence times of superconducting tramsission 
line circuits [H]. We first introduce our model, which 
is based on Ref. [HI As we aim at the realization of a 
pure beam splitter, the Hamiltonian describing our ex- 
perimental system ideally should read as 

n = huj„ {a' a + b'^b) + hg^^ {a'^b + ab^) . (1) 

Here, a) , b\ a, and b are the bosonic creation and an- 
nihilation operators of the two resonators and (7bs is the 
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beam splitter coupling rate. The beam splitter interac- 
tion term ggsici^b + ab^) splits the single resonance sym- 
metrically, resulting in two new normal modes at the an- 
gular frequencies (jJ± =a;o ± (/bs- We can apply Eq. ^ to 
the case of two transmission line resonators coupled in 
a small region around a position where either the mag- 
netic field (current) or the electric field (voltage) has an 
antinode. While this scenario allows one to neglect either 
the capacitive or the inductive coupling channel, it limits 
practical devices to coupling rates smaller than approx- 
imately g / ujaKi2%. In order to achieve higher coupling 
rates, we distribute the coupling over a region spanning 
a significant fraction of the resonator length. As a conse- 
quence of the large coupling strength, the rotating wave 
approximation breaks down, giving rise to a two-mode 
squeezer (TMS) term in the Hamiltonian. Introducing 
the TMS coupling rate ^tms, we find 

TMS (O 

b^ + ab) . 
(2) 

This Hamiltonian describes two coupled harmonic oscil- 
lators of renormalizcd frequency uj, which is split - in gen- 
eral asymmetrically with respect to lOq - into two normal 
modes of frequencies lli± . The detailed definition of w is 
not relevant for this work and can be found in Ref. Uni 
The total coupling rate results from a superposition of 
a capacitive {g^) and an inductive {g\) coupling channel. 
The corresponding two channels couple via 90°-shifted 
single mode fields. Therefore, we find ^bs = (Si + 5c) and 
(7tms = (ft ^ .9c)- The coupling rates gc and g^ depend 
solely on the resonator geometry. For a suitable res- 
onator design, we can achieve gc — gi and thus g-rMs — 0. 
In other words, our distributed coupling approach allows 
for the realization of a pure beam splitter described by 
the Hamiltonian of Eq. ([T]) with large coupling rates (7bs • 

In the remainder of this work, we experimentally vali- 
date the distributed coupling model of Eq. Q by varying 
gc and gi in a controlled way. To this end, we fabricate 
samples containing two coupled microstrip resonators. 
Our design is shown in Fig. [T] For the fabrication of the 
chip shown in Fig. [T](a), we first sputter 100 nm Niobium 
on both sides of a 250 jim thick Si02 (50 nm) coated sil- 
icon wafer. One side is then patterned by optical lithog- 
raphy and reactive ion etching, the other one serves as 
ground plane. Our microstrip waveguides have a width 
of 204 pm to match an impedance of 50 O. As shown in 
Fig.jlja), the two resonators have the same shape. In or- 
der to avoid geometry effects, we shift the position of the 
coupling capacitors defining both ends of the resonators 
rather than redesigning the coupling region. In this 
way, we investigate seven different configurations where 
the resonators are coupled over a length of ^„ = 3 mm at 
different physical coupling positions ii^^ [see Fig. [T|b)]. 
For each two-resonator sample, we fabricate a single res- 
onator sample with the same parameters for comparison. 




FIG. 1. Sample layout, (a) Reworked photograph of two 
coupled resonators on a 10 mm x 6 mm silicon wafer, (b) 
Schematic circuit diagram. The resonators couple over an 
electrical length l^. The coupling region starts at the elec- 
trical length iutt- (c, d) Enlarged view of the region marked 
with the blue (green) box in (a), (e) Sketch of the region 
marked with the purple box in (a). 
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FIG. 2. Transmission spectra of three resonators coupled 
over a length of 3 mm, with 51 = jim (bottom), 300 iim (mid- 
dle), and 600 ym (top). The top and middle curves are shifted 
by 160 dB and 80 dB, respectively. The dashed lines show the 
respective single resonator transmission spectra. 



In our experiments, we measure transmission spectra of 
the fundamental mode of single and coupled microstrip 
resonators at 4.2 K and extract the resonance frequencies 
Wo and a;±, respectively. Typical examples are shown in 
Fig. [2] It can be seen that the two transmission peaks of 
the coupled resonators split asymmetrically with respect 
to the peak of the single resonator. This already indicates 
that the coupling is described by Eq. [2] Furthermore, we 
find that uj„ — 5.44 GHz is independent of the position of 
the coupling capacitors as expected. Hence, also the to- 
tal electrical length f^ot = ttc / {uJo\/(^i) is the same for all 
configurations. With c = 2.99 x 10® m/s and the effective 
dielectric constant tcff — 7.59, we find £tot = 9.963 mm. 
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In order to extract the coupling parameters and g^, 
we make use of the microscopic model |19j leading to 
Eq. ([2]). The input parameters to this model are the 
ratio of the self-inductance (capacitance) per length to 
the mutual inductance (capacitance) per unit length L^^t 
(C„t) and the electrical position of the coupling region. 
The latter cannot be determined directly from the sample 
geometry because the physical length of the transmission 
line differs from the electrical one whenever there is a 
bend in the resonator. Furthermore, the ratios L^^^ and 
C„t depend implicitly on the electrical coupling position. 
Hence, the first step in our analysis is the determination 
of the electrical position of the coupling 



e„ + S£. 



(3) 



Here, as shown in Fig. [T](e), io is the minimum distance 
between the coupling capacitor and the border of the 
coupling region and 6£ accounts for the varying position 
of the coupling capacitor. We obtain 6£ directly from the 
resonator geometry because in good approximation the 
electrical length of a straight segment of the resonator 
equals its physical length. With the definition of Eq. ([S]) 
and the model in Ref. 1191 we can write 



a. 



vcji^o - ^ ) + t^cC^o^+/(^o - 21/^0;^) 
2a;2 



Here, i^l.c = ^oi/Al.c are geometry factors. The 
expressions Al.c represent the overlap integrals 
of the magnetic (electric) field modes. For our 
scenario of homogeneous resonators and funda- 
mental mode coupling, we can set 27rAL c = 

4 T-^tot [sin(27r(C.,-f4)/4ot) - sin(27rA,.,/4„0]- ' In 
order to extract £0 from the measured peak positions Wq 
and Lo±, we first assume that the field in the resonators 
is a TEM-mode and, consequently, L,,,t and C^at are 
independent of ^lett- Subsequent minimization of the 
normalized variance of L^at + Crat for all seven capacitor 
configurations yields £(, = 1.271 mm. Figure [sja) shows 
that indeed for this value of the parameters L^^^ 
and C„t do not deviate more than 3 % from their 
average value. This gives evidence that our model is 
self-consistent. 

In the next step, we use irat, C^^t^ and £„ + 5£ to cal- 
culate (7bs: 5tms, and Co. In Fig. [3]^b), we show g-Bs/^o 
and 5bs/5tms as a function of = Aeft/(4ot - 4)- We 
observe a maximum suppression of the TMS coupling 
rate to 5tms/3bs = 16% and a minimum suppression of 
ffTMs /ffes = 43% while the beam splitter coupling rate 
stays nearly constant at ^bs = (816 ± 7) MHz. An extrap- 
olation of the model prediction suggests that the TMS 
coupling should vanish at the relative coupling position 
£jei = 14%. Nevertheless, the beam splitter coupling rate 
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FIG. 3. Top: The parameters Lrat, Crat, and 4ot divided by 
their respective average value displayed against the relative 
coupling position. Bottom: Coupling ratios as a function of 
the relative coupling position. The solid line is obtained using 
the average values of irat, Crat, and Itot- 



at this position still exceeds 780 MHz. This configuration 
is ideally suited for the realization of a fast beam split- 
ter and can in principle be reached with our geometry. 
Finally, we analyze the potential of our devices for the 
investigation of ultrastrong coupling. In this context, we 
note that the coupling between the two resonators can be 
ultrastrong in the same way as the qubit-resonator cou- 
pling discussed in Ref. [M] For our samples, we achieve a 
maximum TMS coupling rate of 351 MHz for £rei = 27%. 
When moving the coupling region to the center of the res- 
onators, the maximum rate would become 702 MHz and 
(/tms/wq = 12.9%. This implies that the relative coupling 
rate of our device is equally large as in Ref. [TH Since 
our devices do not require nonlinearities and are there- 
fore much easier to fabricate, they provide a promising 
way of studying the dynamics of ultrastrong coupling. 

In summary, we use linear superconducting circuits to 
implement a Hamiltonian with a beam splitter coupling 
strength of more than 800 MHz, where the TMS term is 
suppressed by a factor of six. We demonstrate a tunabil- 
ity of the coupling ratio ^tms/Sbs between 16 % and 43 %. 
An extrapolation of our result shows that an ultrastrong 
coupling scenario as well as a pure beam splitter Hamilto- 
nian can be reached with our sample design. This paves 
the way for studying ultrastrong coupling dynamics and, 
by design of a suitable capacitor configuration [TB] , build- 
ing fast beam splitter circuits for analog quantum com- 
putation and simulation with both standing-wave and 
propagating quantum microwaves. 
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